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CALIBRATING EXTINCTION-FREE STAR FORMATION RATE DIAGNOSTICS WITH 33 GHz FREE-FREE 

EMISSION IN NGC 6946 
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ABSTRACT 

Using free-free emission measured in the Ka-band (26 — 40 GHz) for 10 star-forming regions in the 
nearby galaxy NGC 6946, including its starbursting nucleus, we compare a number of star formation 
rate (SFR) diagnostics that are typically considered to be unaffected by interstellar extinction. These 
diagnostics include non-thermal radio (i.e., 1.4 GHz), total infrared (IR; 8 — 1000 /im), and warm 
dust (i.e., 24 /zm) emission, along with hybrid indicators that attempt to account for obscured and 
unobscured emission from star-forming regions including Ha + 24 /jm and UV + IR measurements. 
The assumption is made that the 33 GHz free-free emission provides the most accurate measure of 
the current SFR. Among the extranuclcar star-forming regions, the 24/itm, Ha + 24 /jm and UV + 
IR SFR calibrations are in good agreement with the 33 GHz free-free SFRs. However, each of the 
SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of 
~2 relative to the 33 GHz free-free SFR. This is more likely the result of excess dust heating through 
an accumulation of non-ionizing stars associated with an extended episode of star formation in the 
nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the 
non-thermal radio continuum yield values which only agree with the 33 GHz free-free SFRs for the 
nucleus, and underestimate the SFRs from the extranuclear star-forming regions by an average factor 
of ^2 and ^4—5 before and after subtracting local background emission, respectively. This result likely 
arises from the CR electrons decaying within the starburst region with negligible escape, whereas the 
transient nature of star formation in the young extranuclear star-forming complexes allows for CR 
electrons to diffuse significantly further than dust heating photons, resulting in an underestimate of 
the true SFR. Finally, we find that the SFRs estimated using the total 33 GHz flux density appear 
to agree well with those from using the free-free emission due to the large thermal fractions present 
at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33 GHz 
observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift 
without the need of ancillary radio data to account for the non-thermal emission. 

Subject headings: cosmic rays - galaxies: individual (NGC 6946) - Hn regions - infrared: general - 
radio continuum: general - stars: formation 
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1. INTRODUCTION 

Accurate estimates of star formation rates (SFRs) are 
crucial in nearly all fields of extragalactic astronomy. Be- 
sides requiring a high level of accuracy, the ease with 
which such estimates can be made plays a critical role 
in the overall utility of a SFR diagnostic given that, 
for studying galaxy evolution, observations need to be 
made over large areas of sky and be deep enough to de- 
tect galaxies at high redshifts. Presently, a wide variety 
of SF R diagnostics a re commonly used in the literature 
(e.g. JKeim icutt 1998|), each of which have their strengths 
and shortcomings. 

Rest-frame UV (1250 - 2500A) and optical (e.g., Ha) 
observations provide a measure of the radiant energy re- 
leased by newly formed massive stars. While directly 
related to the massive star formation process, such ob- 
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servations are hampered by interstellar extinction that 
varies enormously within and among galaxies. This is 
especially problematic for trying to characterize the evo- 
lution of the SFR density with redshift as the galax- 
ies that dominate the stellar mass assembly between 
1 ^ z 3 appear to be dusty starbursts with in- 
frared (IR; 8 — 1000 /im) luminosities ranging between 
10 10 Lp, < £ T r < 10 13 (e.g.. ICharv fc Elbad [200ll: 
Takeuchi etlhll2005tlBouwens et al.ll2009t iMurphv et alJ 



20111) . In these IR-bright starbursts, mor e than >90% 



of the UV emission can be obscured (e.g., iHowell et al.l 
120101 IBuat eTaII[20T0l: IMurphv et al.ll201lh . 

Rest-frame IR measurements, particularly those in the 
far-infrared (FIR; 42.5 — 122.5 /im), do not suffer sig- 
nificantly from extinction as they probe the re-radiated 
UV/optical emission from heated dust grains surround- 
ing massive star-forming regions. While an extinction- 
free measure, the interpretation of IR emission alone 
is also complex. Variations in the dust composition, 
content, and distribution along the line of sight will 
affect the fraction of UV photons absorbed, while a 
portion of the IR emission will arise from dust heated 
by older stars (i.e., the cirrus component | Heloul 



Lonsdale Persson fe Helou 1987; Hirashita et al.l 



1986 



2003 



Bendo et al.l l2010). 



Radio continuum emission is also widely used as a 
tracer of star formation in galaxies both at low and 
high redshift as the result of the tight empirical corre- 
lation between their optically thin synchr otron emission 
at 1. 4 GHz and their FI R radiation (e.g.. Ide Jong et aLI 
119851: iHelou et al.l 11985ft . While this relation (i.e., the 
FIR-radio correlation) is certainly rooted in common de- 
pendencies on massive star formation, it is unclear how 
presumably unrelated physical processes affecting the 
propagation of cosmic-ray (CR) electrons and the heat- 
ing of dust grains work together to yield a nearly ubiq- 
uitous corr elation over near ly 5 orders of magnitude in 
luminosity (|Yun et al.l [200 1). Spitzer revealed local cor- 
relations in the spatial distributions of 70 /im and non- 
thermal radio emission in the disks of galaxies that reflect 
an age effect indicating the complexities of using non- 
thermal radio continuum em ission as a SFR diagnostic 
(jMurphv et al.ll2006bl 120081 ). Additional complications 
for using the non-thermal radio emission as a SFR diag- 
nostic are also evident given that the correlation between 
the IR and radio emission is strikingly different for spiral 
arms and interarm regions (|Dumas et al.|[20TlT) . 

Higher frequency radio continuum emission, which be- 
comes dominated by thermal (free-free) radiation once 
beyond ~30 GHz for globally i ntegrated mea surements 
of star-forming galaxies (e.g., iCondonl Il992| ). is both 
largely extinction free and can be directly related to 
the ionizing photon rate arising from newly formed mas- 
sive stars. Evidence for this has been indicated by 
detail ed studies of star-forming regions in the Galaxy 
(e.g., IMezger fc Henderson! [l967Tl . nearby dwarf irreg- 
ulars (e.g. iKlein fe Graevd 119861). the nuc lei of nor- 
mal galaxie s (e.g.. | Turner fc Hoi 119831 I1994D a nd star- 
bursts (e.g.. lKlein et al.ll988HTurner fc Holll985l) . as well 
as high resolution investigations of super star clusters 
within nearby blue compact dwar fs (e.g., iTurner et al.1 
1998; IKobulnickv fc Johnsonl[T999t) . 

While typically considered an ideal measure for the 
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Fig. 1. — The 24 fim map of NGC 6946 showing the location of 
each star-forming region in the present investigation. The diameter 
of each circle corresponds to 25", the FWHM of the GBT beam 
at 33 GHz, which sets the resolution of this study. This diameter 
projects to a linear scale of 0.8 kpc at the distance of NGC 6946. 



current star formation activity in galaxies, the presence 
of an 'anomalous' dust emission component in excess of 
the free- free emission between ~10 and 90 GHz, generally 
attributed to rapidly rotating ultrasmall grains with a 
non-zero electric moment (e.g.. | Draine fc Lazarianlll998l : 
iPlanck Collaboration et all 1201 lh . may complicate this 
picture. F or a single outer-disk star-forming region in 
NGC 6946. IMurphv et al.l (l2l)10l) detected excess 33 GHz 
emission relative to what is expected given existing lower 
frequency radio data. This result is interpreted as the 
first likely detection of so-called anomalous dust emission 
outside of the Milky Way. Given that the excess was only 
detected for a single region, this emission component may 
be negligible for globally integrated measurements. 

In this paper we quantitatively compare SFR estimates 
using a number of "extinction-free" diagnostics for 10 
star-forming regions within the n earby galaxy NGC 6946 
presented in IMurphv et al. (2010). This is done to assess 
their reliability against SFRs measured via free-free emis- 
sion at 33 GHz. The paper is organized as follows: In §2 
we introduce the data sets used for our SFR comparisons. 
Then, in §3, we calibrate each of the observed quantities 
for the same initial mass function (IMF). Comparisons 
between each of the SFR diagnostics are made in §4 and 
discussed in §5. Our main conclusions are summarized 
in §6. 

2. DATA AND PHOTOMETRY 

NGC 6946 is a nea rby (d 6.8 Mpc; 

iKarachentsev et ail I2000D late-type s piral galaxy 
harboring a mild s tarburst n ucleus ([Turner fe Hoi 
[1981 IBall etaLl fl985h . Recently. IMurphv et al.1 (2010) 
targeted 10 star-forming regions within NGC 6946, 
including the nucleus, using the 100 m Robert C. 
Byrd Green Bank Telescope (GBT) in the Ka-band 
(26 — 40 GHz). These regions were selected because they 
have existing and forthcoming mid- and far-infrared 
spectroscopic data collected as part o f the Spitzer In- 
frared Nearby Galaxies Survey (SINGS: iKennicutt et al.1 
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TABLE 1 

Source Positions and Radio Photometry 



R.A. Decl. 533 GHz Sg,s GHz S4.9 GHz Si. 7 GHz Si. 5 GHz Si. 4 GHz 

(J2000) (J2000) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) 



Nucleus 


20 34 52.34 


+60 9 14.2 


15.7+0.79 


31.8+1.59 


44.8+2.24 


77.6+3.88 


83.2+1.59 


85.5+1.59 


Enuc. 


1 


20 35 16.65 


+60 11 1.1 


1.0t±0.05 


0.9+0.12 


0.9+0.09 


1.1+0.15 


1.1+0.12 


1.3+0.12 


Enuc. 


2 


20 35 25.49 


+60 10 1.8 


2.4+0.16 


2.6+0.17 


3.1+0.17 


4.3+0.26 


4.5+0.17 


5.0+0.17 


Enuc. 


3 


20 34 51.89 


+60 12 44.8 


1.0+0.07 


0.9+0.12 


1.1+0.10 


1.6+0.16 


1.7+0.12 


1.8+0.12 


Enuc. 


4 


20 34 19.17 


+60 10 8.7 


2.9+0.15 


1.4+0.13 


1.8+0.12 


2.0+0.17 


2.1+0.13 


2.3+0.13 


Enuc. 


5 


20 34 39.27 


+60 4 55.1 


0.6t±0.05 


0.6+0.11 


0.6+0.09 


1.0+0.15 


1.0+0.11 


1.2+0.11 


Enuc. 


6 


20 35 6.09 


+60 11 0.6 


2.8+0.18 


3.6+0.21 


4.3+0.23 


6.7+0.37 


7.2+0.21 


7.6+0.21 


Enuc. 


7 


20 35 11.21 


+60 8 59.7 


3.1+0.20 


2.9+0.18 


3.6+0.20 


4.8+0.28 


5.1+0.18 


5.3+0.18 


Enuc. 


8 


20 34 32.52 


+60 10 22.0 


2.2t±0.14 


2.0+0.15 


2.4+0.14 


2.9+0.20 


3.1+0.15 


3.5+0.15 


Enuc. 


9 


20 35 12.71 


+60 8 52.8 


2.4+0.16 


2.5+0.17 


2.9+0.17 


4.2+0.25 


4.4+0.17 


4.9+0.17 



Note. — ' Correction ('boosting') factors arising from an over subtraction of signal by reference beams have been 
applied to 33 GHz flux densities for sources in which the correction was larger than 15%. These scale factors are 2.55, 
1.17, and 1.25 for Enuc. 1, 5, and 8, respectively. 



TABLE 2 
Infrared Photometry 



ID 




/„ (8 M m)t 


U (24 M m) 


U (70 Atm) 


/„ (100 Mm) 


f v (160 A*m) 


/„ (250 Mm) 


U (350 A*m) 


U (850 Aim) 






(mjy) 


(mjy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(mjy) 


Nucleus 


1506.7+ 75.4 


4308.9+ 215.4 


26.11+ 2.61 


32.09+ 3.21 


23.67+ 4.74 


6.80+ 1.09 


3.12+ 0.50 


401.8+ 62.8 


Enuc. 


1 


33.6± 2.8 


88.0+ 4.5 


0.85+ 0.13 


1.23+ 0.18 


1.02+ 0.23 


0.33+ 0.06 


0.15+ 0.04 


10.8+ 17.9 


Enuc. 


2 


53.0+ 3.5 


175.1+ 8.8 


1.28+ 0.16 


1.94+ 0.24 


1.65+ 0.34 


0.48+ 0.08 


0.22+ 0.05 




Enuc. 


3 


40.6+ 3.1 


66.4+ 3.5 


1.00+ 0.14 


1.54+ 0.21 


1.30+ 0.28 


0.36+ 0.07 


0.20+ 0.05 




Enuc. 


4 


51.9+ 3.5 


83.1+ 4.3 


1.00+ 0.14 


1.65+ 0.21 


1.49+ 0.32 


0.35+ 0.07 


0.17+ 0.04 


15.2+ 17.9 


Enuc. 


5 


12.5+ 2.4 


24.7+ 1.6 


0.40+ 0.11 


0.57+ 0.15 


0.60+ 0.16 


0.16+ 0.04 


0.11+ 0.04 




Enuc. 


6 


192.3+ 9.9 


352.2+ 17.6 


2.62+ 0.28 


3.77+ 0.40 


3.02+ 0.61 


0.69+ 0.12 


0.37+ 0.07 


54.7+ 19.6 


Enuc. 


7 


186.6+ 9.6 


366.0+ 18.3 


3.21+ 0.34 


4.37+ 0.46 


3.65+ 0.74 


1.21+ 0.20 


0.57+ 0.10 


60.5+ 20.0 


Enuc. 


8 


114.7+ 6.2 


178.5+ 9.0 


1.55+ 0.19 


2.27+ 0.26 


1.92+ 0.40 


0.56+ 0.10 


0.26+ 0.05 


39.6+ 18.8 


Enuc. 


9 


170.3+ 8.8 


301.2+ 15.1 


2.25+ 0.25 


3.31+ 0.36 


3.03+ 0.61 


0.95+ 0.16 


0.46+ 0.08 


58.8+ 19.9 


Note. 




^ Corrected for stellar emission. 


Observations 


were made using 


the following 


observatories and 


! instruments: 


Spitzer-IRAC 



(8 Aim), Spitzer-MlPS (24 Aim), HerscheLPACS (70, 100, 160 Aim), HerscheLSPIRE (250, 350 A«m), and SCUBA (850 A»m). 



TABLE 3 
Ha and UV Photometry 



ID 




f 1 " 


/„ (1528 A) 


/„ (2271 A) 


E(B- V)* 






(10 -12 ergs -1 cm -2 ) 


(mjy) 


(mjy) 


(mag) 


Nucleus 


1.09+ 0.22 


0.07+ 0.01 


0.42+ 0.06 


0.342 


Enuc. 


1 


0.47+ 0.09 


1.25+ 0.19 


2.49+ 0.37 


0.345 


Enuc. 


2 


2.01+ 0.40 


5.81+ 0.87 


9.58+ 1.44 


0.346 


Enuc. 


3 


0.80+ 0.16 


2.52+ 0.38 


4.44+ 0.67 


0.343 


Enuc. 


1 


0.61+ 0.12 


1.65+ 0.25 


3.30+ 0.50 


0.343 


Enuc. 


5 


0.45+ 0.09 


1.04+ 0.16 


2.06+ 0.31 


0.335 


Enuc. 


6 


1.09+ 0.22 


1.80+ 0.27 


3.75+ 0.56 


0.344 


Enuc. 


7 


1.19+ 0.24 


1.93+ 0.29 


4.07+ 0.61 


0.342 


Enuc. 


8 


0.36+ 0.07 


0.23+ 0.04 


0.67+ 0.10 


0.342 


Enuc. 


9 


0.76+ 0.15 


1.16+ 0.17 


2.44+ 0.37 


0.342 



Not e. — t Ha imaging obtained wi th the KPNO 0.9 m by IFerguson et all 
dl998D . * Galactic extinction taken from (Schlcgcl ct al. 1998) used to correct the 
Ha and GALEX FUV and NUV flux densities assuming A v /E ( B - V) = 3.1 a nd 
the modeled extinction curves o f lWeingartner fe D rainc (2001); Drainc 120031) . 



2003) and the project, Key Insights on Nearby Galaxies: 
a Far-Infrared Survey with Herschel (KINGFISH; PI. 
R.C. Kennicutt). The position of each source is given in 
Table [1] In Figure [T] we show the location of each region 
on the 24 /j,m image of NGC 6946; the corresponding cir- 
cle diameters are 25", which matches the FWHM of our 
lowest resolution data for the present multiwavelength 
SFR comparison. This projects to a physical scale of 
« 0.8 kpc at the distance of NGC 6946. 



2.1. Radio Data 

Radio data at 1.4, 1.5, 1.7, 4.9, 8.5, an d 33 GHz used 
here a re the same as the data presented in lMurphv et al.l 
(2010), where a more detailed description of the obser- 
vations and spectral decomposition of the thermal and 
non-thermal components can be found. The 1.4 GHz 
radio map (14" x 12"5 beam) comes from the Wester- 
bork Synthesis Ra dio Telescope (WSRT)-SINGS survey 
(jBraun et al.ll2007t) . Imaging at 1.5, 1.7, 4.9, and 8.5 GHz 
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(15" x 15" beam) all come from|BecJj (|2007l ). We note 
that the 4.9 and 8.5 GHz radio data included single-dish 
measurements for short-spacing corrections. 

Observations at 33 GHz were taken on 2009 March 21 
using the Caltech Continuum Backend (CCB) on the 
GBT. The CCB simultaneously measures the entire Ka 
bandwidth over 4 equally spaced frequency channels. 
The beamswitched signal is synchronously read out and 
demodulated to remove atmospheric fluctuation and/or 
gain variations. Reference beams are measured by nod- 
ding 1 '3 away from the source; details abo ut the reference 
beam locations can be found in Figure 1 of lMurphv et al.1 
(2010). The average FWHM of the GBT beam in the 
Ka-band was «25" among our sets of observations. A 
detailed description on the performance of the CCB re- 
ceiver, the data reduction pipe line, and error estimates 
are given m IMason et all (|2009| ). 

2.2. UV, Optical, and IR Data 

GALEX far-UV (FUV; 1528 A) and near-UV (NUV; 
2271 A) data of NGC 6946 were taken from the GALEX 
archive and will be included in the GALEX Large Galaxy 
Atlas (M. Seibert et al. 2011, in preparation). The 
calibration uncertainty for these data is ~15% in both 
bands. The Ha image of NGC 6946, obtained using the 
Tek 2048 x 2048 CCD on the KPNO 0. 9 m (23' x 23' 
FOV and 0'.'69 pixel -1 ), was taken from Ferguson et al.1 
(1998), and has a calibration uncertainty of «20%. The 
8, 24, and 850 /im data used here were included in the 
SINGS fifth data release, 18 and have calibration uncer- 
tainties of 5, 5, and 15%, respectively. Regions 2, 3, and 
5 were not covered by the SCUBA 850 /im maps, while 
regions 4 and 8 were near the edges of the maps. We note 
that the 8 /im map was corre cted for stellar light using a 
scaled 3.6 /im imag e following iHelou et all (pool . We do 
not use the Spitzer 70 /im or the SCUBA 450 /im data, 
also included in the SINGS fifth data release, in lieu of 
new Herschel data taken as pa rt of KINGFISH. 

We make use of the Herschel (|Pilbratt et al.l l2010) data 
obtained by the KINGFISH project having resolutions 
better than the 25" beam of our GBT data which in- 
clude PACS 70, 100, and 160 /zm imaging, along with 
SPIRE 250 and 350 /im imaging. Both PACS and SPIRE 
imaging were reduced using t he Hersch el Interactive Pro- 
cessing Environment (HIPE; Ott 2010) software package, 
version 3.0.1510. The current calibration has uncertain- 
ties of -10, 10, an d 20% for the 70, 10 and 160 /xm 
bands, respectively (jPoglitsch c t al.1 12010ft . The uncer- 
tainty in the SPIRE 250 and 350 /im data is taken to 
be 15%, which is the quadrature sum of the 15% abso- 
lute calibration uncertainty and a 1% uncertainty in the 
size of the beam; following the current recommendation 
by the SPIRE Observer's Manual, we assume beam sizes 
at 250 and 350 /xm of 426 and 771 arcsec 2 , respectively. 
Calibration correction factors of 1.02 and 1.05 were ap- 
plied to the 250, 350 /im data as suggested in the SPIRE 
Scan-Map AOT Observer's manual. 

2.3. Photometry 

The photometry was carried out on all UV, optical, in- 
frared, and radio maps after first cropping each image to 

18 SINGS data products can be found at 
|http://irsa.ipac.caltech.edu /data/SPITZER/SINGS / 



a common field-of-view and regridding to a common pixel 
scale. To accurately match the photometry of our images 
to the GBT measurements, maps were convolved to the 
resolution of the GBT beam in the Ka- band following 
the im age registration description given in lMurphv et al.1 
(2010). Registration of the Herschel data was carried out 
in a similar way as the registration of the Spitzer data. 
Each of the Herschel images were first convolved to the 
resolution of the 350 / im data using custo m smoothing 
kernels as described in lGordon et al.l ()2008l ) and updated 
by those authors for Herschel, CLEANcd to suppress 
side-lobe structure in the PSF, and then restored using a 
Gaussian beam having a FWHM of 25". Flux densities 
at each wavelength were measured by taking the surface 
brightness at the location of each GBT pointing and mul- 
tiplying by the effective area of the beam. In the case 
of the UV and Ha photometr y, we correct e a ch re gion 
for Milky Way extinction using ISchlegel et all (fl998h as- 
suming Ay/E{B — V) = 3.1 and the m o deled extincti on 
curves o f iWeingartner fc Drainel (|200lh ; iDralnel pOOl . 

Unlike IMurphv et al.l (12010ft . we subtract a local back- 
ground estimate from our photometry within the vicinity 
of each star-forming region to remove any diffuse emis- 
sion component that is most likely unassociated with on- 
going star formation. Results using the photometry with- 
out subtracting local backgrounds are discussed in ^15.41 
and in the Appendix. The local background estimates 
were measured by placing four 25" diameter apertures 
at a distance of 1.5 times the beam FWHM (i.e., 37'.'5) 
away from the center of the source positions in each of the 
four Cardinal directions. The median surface brightness 
among the pixels within these four apertures (typically 
100 pixels in total) was then multiplied by the effective 
area of the beam to get an estimate of the local diffuse 
background emission. While such a local background 
subtraction for our 33 GHz GBT single-pointing obser- 
vations is not possible, we note that those observations 
were taken using 1.'3 nods, thus reference beams were 
typically measured on diffuse parts of the galaxy and 
likely provide a reasonable measure of the local back- 
ground. Furthermore, at 33 GHz, we do not expect there 
to be a significant underlying diffuse component if the 
emission is in fact dominated by free-free emission. 

Photometric uncertainties were conservatively esti- 
mated by taking the quadrature sum of the calibration 
and background uncertainties, along with the RMS noise 
for each image. The results from our radio, infrared, and 
UV/Ha photometry, along with \-a uncertainties, are 
given in Tables [TJ [21 and |31 respectively. The local back- 
ground values, along with the average fractional contri- 
bution of the background emission at each waveband and 
region, are given in the Appendix. For the cases in which 
an over subtraction of signal in the 33 GHz observations 
has been estimated to be at the >15% level due to refer- 
ence beams landing on bright regions in the galaxy disk, 
we apply correcti on ('boosting') facto rs to the 33 GHz 
flux densities fsee IMurphv et aLll2010D . These scale fac- 
tors were derived based on the 8.5 GHz and 24 /im im- 
ages, and are approximately 2.55, 1.17, and 1.25 for 
Enuc. 1, 5, and 8, respectively. We note th at these values 
are < 2% smaller than those reported by IMurphv et al.l 
( 2010) due to including local background subtractions. 
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2.4. Modeling the Infrared-to- Radio Spectra 

We fit the radio and infrared data independently given 
that the physical processes producing these emissions are 
distinct. The infrared (i.e., 8 — 850 /im) photometry are 
fit by the spe c tral e nergy distribution (SED) models of 
iDale fc Heloul ([20021) : associated total infrared (IR: 8 - 
1000 /im) luminosities from integrating the best-fit SEDs 
for each region are given in Table 2] Fitting errors in 
the IR luminosity determinations were estimated by a 
standard Monte Carlo approach using the photometric 
uncertainties of the input flux densities. 

The radio data are fit by varying a combination of 
thermal (free-free) and non-thermal (synchrotron) emis- 
sion components which scale as Sj oc i/~ aT and S^ T oc 
^° NT , respectively, where ax = 0.1 and «nt (see Ta- 
ble |4] for estimated values) are the thermal and non- 
thermal spectral indices, respectively. The only excep- 
tion to this is for extranuclear region 4 where the ther- 
mal and non-thermal emission were fit only to the radio 
data at v <C 10 GHz and an additional component (i.e., 
spinning dust emission; lAli-Haimoud et all 120 09) was 
used to fit the excess 33 GHz emission (see lMurphv et al.l 
120101 for details) . Recent follow-up observations of this 
region at 16.5 GHz using the Arcminutc Microkelvin Im- 
ager have confirmed a rising spectrum between 8.5 and 
15 GHz ([Scaife et al.l [2010). The excess emission ac- 
counts for «56% of the total 33 GHz flux density from 
extranuclear region 4 and has been subtracted out for 
the present analysis. We note that the quantitative re- 
sults prese nted here are s i gnific antly different from what 
is given in iMurphv et al.l (|2010t) as those authors did not 
subtract out a local background estimate and included a 
spinning dust component for star-forming regions in ad- 
dition to extranuclear region 4. Results associated with 
fitting the photometry without subtracting a local back- 
ground can be found in §5.41 and the Appendix. 

Before varying the free-free and non-thermal compo- 
nents to fit the observed radio spectra, we first estimate 
the non-thermal spectral index for each source. The 
non-thermal emission is proportional to the ratio of en- 
ergy losses of CR electrons from synchrotron emission 
to total energy losses. In our model, we account for 
energy losses to CR electrons arising from synchrotron 
radiation, inverse Compton (IC) scattering, ionization, 
bremsstrahlung, and escape through an empirical pre- 
scription. Additional losses through advection from a 
galactic scale wind are ignored. A much more detailed 
description on the con struction of th e modeled radio 
spectra can be found in iMurpriyl (|2009| ). 

Synchrotron losses are proportional to Ub = B 2 /(8ir), 
the magnetic field energy density, while IC losses are pro- 
portional to [/ ra d, the radiation field energy density. We 
use magnetic field strengths estimated from the 1.4 GHz 
radio continuum surface brightnesses f ollowing the re- 
vised minimum energy calculation of iBeck fc Krausd 
(|2005l see Table [4]). For this calculation we assume a 
proton-to-electron number density ratio of K w 100 ±50 
and a pathlength of I w 1 ± 0.5 kpc. The estimated 
-Bmin values range between w 10 — 16 /xG among the 
extranuclear star-forming regions and it is w 30 /iG 
for the nucleus. Assuming that the bolometric sur- 
face brightness of these dusty star-forming regions is 
well approximated by the IR surface brightness, we es- 



timate ?7 la d ~ (27t/c)/ir. These values range between 
w 0.4 — 4.2 x 10 -12 ergs cm -3 among the extranuclear 
star- forming regions, and it is « 41 x 10~ 12 ergs cm" 3 
in the nucleus. We assume that the density of the in- 
terstellar medium (ISM) is 7Jism = 0.1 and 10 cm~ 3 for 
the extranuclear regions and nucleus, respectively. Tak- 
ing these values, we make a rough estimate for the shape 
(i.e., spectral index) of the non-thermal radio spectrum. 
Once the non-thermal spectral index is fixed, the thermal 
and non-thermal components are varied to fit the radio 
observations. 

The radio to infrared spectra of each star-forming re- 
gion are plotted in Figure [21 along with the components 
of the fits. In Table [4] we list the observed (a obs ) and 
(estimated) non-thermal radio spectral indices measured 
between 1.4 and 8.5 GHz, along with the associated ther- 
mal radio fractions at 33 GHz from our radio spectral fit- 
ting for each star-forming region. The observed indices 
are quite flat, being «0.2, on average, consistent with 
young Hn regions, while the average non-thermal spec- 
tral index is «0.8. The average 33 GHz thermal fraction 
is found to be ~87% among these star-forming regions. 

3. SFR CALIBRATIONS 

To ensure that each SFR diagnostic can be compared 
fairly, w e derive correspondin g calibrations using Star- 
burst99 dLeitherer et al.lll999D for a common IMF. We 
choose a Kroupa ([Kroupall2001| ) IMF, having a slope of 
— 1.3 for stellar masses between 0.1 — 0.5 Mq and —2.3 
for stellar masses ranging between 0.5 — 100 Mq . Thes e 
calibrations update those presented in lKennicuttl (|1998|h 
Only SFR diagnostics that are considered to be reliable 
independent of extinction are discussed. 

3.1. UV, Optical, and IR SFRs 

Assuming a solar metallicity and continuous star for- 
mation (i.e., a fixed SFR), Starburst99 stellar population 
models yield the following relation between the SFR and 
production rate of ionizing photons, Q(i?°), at an age of 
-100 Myr: 



SFR 



M Q yr 



7.29 x 10 



-54 



Q(H°V 



(i) 



For reference, the coefficie nt relating t he SF R and ion- 
izing photon rate given in iKennicuttl (11998|) is a factor 
of pal . 5 times larger than that above (see ICalzetti et all 
120071 for a discussion on the origin of this difference). 
We choose a constant SFR over a timescale of —100 Myr 
for all estimators as we are sampling star-forming com- 
plexes on the scales of «0.8 kpc that likely include a range 
of ages among individual star formation sites. Conse- 
quently, below certain spatial scales (i.e., < a few hun- 
dred pc; the size of individual giant Hn regions) these 
calibrations will likely break down as the assumption of 
spatial averaging is no longer valid. 

The ionizing photon rate can of course be expressed as 
a (extinction corrected) H recombination line flux, such 
that for Case B recombination, and assuming an elec- 
tron temperature T e = 10 K, the Ha recombination line 
strength is related to the SFR by 



' SFRho 

,Mq yr- 



= 5.37 x 10~ 42 



L 



'He 



erg s 



-l 



(2) 
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tween SFR and FUV luminosity, 



10° 10' 10 s 10 3 10" 10° 10' 10 2 10 3 10" 
v (GHz) 

Fig. 2. — The radio to infrared spectra of each star- forming 
region. The radio and infrared photometry are plotted as filled 
circles; error bars are shown, but are usually smaller than the 
size of the plotting symbol. F its to the dat a using a combina- 
tion of mo deled radio spectra {M urphy 2009) and infrared SED 
libraries of Dale & Hclou (2002) are shown as solid lines. The 
infrared dust emission is shown by a short-dashed line. The 
radio spectrum is modeled by 2 components, non-thermal syn- 
chrotron emission (long-dashed lines) and thermal (free-free) emis- 
sion (triple-dot dashed line) except for the case of extranuclear 
regio n 4, where an additional component (i.e., spinning dust emis- 
sion; [Ali-Hai'moud ct al. 2003, dotted line) was ne eded to fit the 
excess emission at 33 GHz (sec Murphy ct al. 2010, for details). 



The above equation indicates that the SFR is directly 
proportional to the Ha line luminosity, assuming that a 
constant fraction of the ionized H atoms will emit an Ha 
photon as they recombine, and that the extinction cor- 
rection is accurate. However, if a significant fraction of 
ionizing photons are absorbed by dust, the above equa- 
tion will underestimate the SFR, and the = sign should 
be replaced by > . Since the ionizing flux comes from very 
massive stars with lifetimes <10 Myr, we note that the 
coefficients in Equations [T] and [2] are nearly independent 
of starburst age under the assumption of continuous star 
formation so long at it is >10 Myr. Accordingly, such 
measurements sample the current (i.e., ~10 Myr) star 
formation activity. 

The integrated UV spectrum is dominated by 
young stars, ma king it a sensit i ve probe of recent 
(~10 - 100 Myr; IKejmicutt] Gill ICalzetti et afl [20051: 
ISalim et al.l 120071) star formation activity. We convolve 
the output Starburst99 spectrum with the G ALEX FW 
transmission curve to obtain the following conversion be- 



/ SFR 



FUV 



VM Q yr- 



= 4.42 x 10" 



-^FUV 

erg s _1 



(3) 



To derive a calibration for the total infrared (IR; 
8 — 1000 /xm), we make the assumption that the entire 
Balmer continuum is absorbed and re-radiated by dust 
and that the dust emission is optically thin. Integrating 
the output Starburst99 spectrum over this wavelength 
range (i.e., 912 < A < 3646 A) results in the following 
(predicted) relation between the IR emission and current 
SFR: 



SFR 



IR 



Mq yr 



= 3.88 x 10" 



L 



in 



erg s 



(4) 



Since a non-negligible fraction of the far-infrared con- 
tinuu m may be heated by an older stellar population 
(e.g., ISau vagc fc Thuanl 1992; Walterbos fc Greenawaltl 
ll996tlBendo et al.ll2010fl . it is likely more physical to only 
relate warmer du st emission to the current SFR (e.g., 
IHelou et al.1 12004). Thus, we also include a monochro- 
matic 24 /im based SFR diagnostic, for which there are 



a number of calibrations in the literature (e.g-JWu et al 
2005; iPerez-Gonzalez et al.l 120061: fA lonso-H errero et al 
20061: ICalzetti et al.ll2007t iRelano et al.ll2007t iZhu et al 



2008; Ric ke et al.l 12009). A detailed com parison among 



-A 

each o f these relations can be found in ICalzetti et al.l 
(2010). For simplicity, we only cho ose a single calibra- 
tion to use in our comparison (i.e., Rel ano et al.l l200~7l 
calibrated for a Kroupa IMF) , such that 



= 5.58 x 10- 36 

M Q yr-i J 



vL„(2A /im) 



-, 0.826 



erg s 



(5) 



as this relation has been calibrated over the largest 
range in 24 /im luminosity (i.e. 1 x 10 38 erg s 



-i < 



vL v {2A urn) < 3 x 10 44 erg s" 1 ). 

Given that not all of the UV/optical photons will be 
absorbed and re-radiated by dust, a series of new em- 
pirical calibrations based on the linear combination of 
observed 24 /im (obscured star formation) and Ha (unob- 
scure d star formation) luminosities have been developed 
(e.g.. ICalzetti et al.l 120071 : iKennicutt et al l 120071 120091: 
IZhu et al.ll2008l ). Based on the analysis of ICalzetti et al.l 
(|2007f ). we adopt the following relation for our sample of 
star-forming regions to obtain extinction corrected Ha 
line luminosities: 



j corr 

erg s _] 



r obs 



0.031^(24 /im) 



erg s 



(6) 



While this relation has been calibrated for Hn regions, 
we note that it has also been found to be appropriate 
for entire galaxies having 4 x 10 42 < i/L„(24: /im) < 5 x 
10 43 erg s _1 . For galaxies having vL v {2A /im) < 4 x 
10 42 erg s _1 , this relation is optimized by rep lacing the 
coefficient 0.031 by 0.020 (jCalzetti et al.H201fjD . which is 
th e calibration derived fo r normal star-forming galaxies 
by IKennicutt et al~l ([2009D . Combining Equation [5] with 
Equation [5] results in the following SFR relation, 



SFR 



M Q yr 



mix * =5.37xl0" 42 



t obs 



0.031iAL„(24 /tm) 



erg s 



(7) 
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Similarly, by combining the SFR estimates from the 
total IR and (observed) UV emission, one can account 
for the obscured and unobscured emission contributing 
to the total (bolometric) SFR, which we define as: 



SFRt, 



SFRfuv + SFRir. 



(8) 



Using the above calibrations, this relation can be ex- 
pressed as a linear combination of the UV and total IR 
emission such that, 



SFRt °Q = 4.42 x lCr 44 
M Q yr"V 



L 



FUV 



0.88L m 



erg s 



(9) 



This diagnostic is often used to c haracterize SFRs from 



20061: iBuat et al.l 12007 


, 120111) and high- 


Iglcsias-Paramo ct al. 


20041: lElbaz et al.l 


Daddi ct al. 2007: Rcddv ct al. 2010) studies. 



2007; 



as already mentioned above, the presented calibrations 
will depend on each complex including a range of ages 
among individual star formation sites. Therefore, cali- 
brations such as these will have increasing uncertainties 
when applied to a single star-forming region depending 
on how well the IR and UV are measuring emission 
from recent star formation, not to mention the fact 
that the IR and UV are sensitive to star formation on 
different timescales (i.e., ^10 Myr versus ~10— 100 Myr, 
respectively). It is also worth noting that the derived 
coefficient for scaling the IR luminosity given above 
is larger than that from empirical studies of starburst 
galaxies which report a coefficien t closer to ~0.6 (e.g., 
IMeurer et alllllM ICalzettil[200l . 



3.2. Radio SFRs 

Radio continuum emission from galaxies is gener- 
ally composed of two optically-thin components: non- 
thermal synchrotron emission associated with CR elec- 
trons accelerated in a galaxy's magnetic field, and ther- 
mal bremsstrahlung (free-free) emission around massive 
star-forming regions. The origin of each of these compo- 
nents lies in the process of massive star formation. 

At high radio frequencies, where r <C 1, the ionizing 
photon rate is directly proportional to the thermal spec- 
tral luminosity, Lj, , varying only weakly with electron 
temperature T e |Rubin 1961), such that 



Q(H°) 



6.3 x 10 25 



T 



10 4 K 



-0.45 



GHz 



L 



erg s 1 Hz 



(10) 



As with the H recombination line fluxes, it is again worth 
noting that Q(H°), and consequently the SFR, may in 
fact be underestimated by the free-free emission if a sig- 
nificant fraction of ionizing photons are absorbed by dust; 
in this case the = sign in the above equation should be 
replaced by >. By combining Equations [Tl and 1101 one 
can derive a relation between the SFR and thermal radio 



emission: 
SFR^ 

M q yr" 



= 4.6 x 10~ 28 



10 4 K 



-0.45 



Vghz/ 



v erg s 



-i Hz- 1 
(11) 

Using a Kroupa IMF that extends down to 0.1 Mq results 
in a coefficient that is a factor of ~2.5 times larger than 
that found in the relation between t he S FR and free- 
free radio emission given in ICondonl (U992T). an d ~25% 
smaller than that given in iSchmitt et aTTpOOl ). 

At lower radio frequencies, which are typically domi- 
nated by non-thermal synchrotron emission, calibrations 
between the supernova rate, and thus the SFR, have been 
developed. From the output of Starburst99, which as- 
sumed a supernova cut-off mass of 8M Q , we find that 
the total core-collapse supernova rate, Asn, is related to 
the star formation rate by, 



SFR 



M Q yr- 



86.3 



yr" 



(12) 



Work comparing the non-thermal spectral luminosity 
with the supernova rate in the Galaxy has yielded an 
empirical calibration such that 



SN 



GHz 



= 1.3 x 10 30 , — 
v erg s 1 Hz 1 J I yr 

^_ , , , ' (13) 

(jTammann 1982; Condon fc YirJ I1990T ). By combining 
Equations [12] and 1131 we can express the star forma- 
tion rate as a function of the non-thermal radio emission 
where 



SFR 



NT 



Mq yr" 



= 6.64xl0 -29 



/ v 



GHz 



L 



NT 



erg s 



Hz- 



(14) 

We note that using the Kroupa IMF that extends down 
to 0.1 Mq results in a coefficient that is a factor of ~3.5 
times larger than that used in the expres sion relating th e 
SFR and non-thermal radio emission bv ICondonl (fl9 92). 
and_~50% smaller than that given in ISchmitt et al.l 
<f2006h . 

Since the observed radio continuum emission is com- 
prised of both free-free and synchrotron emission, we 
can combine Equations [11] and Q3] to construct a single 
expression for the SFR from the total radio continuum 
emission at a given frequency such that: 



SFR„ 

Mq yr- 



10 



-27 



2.18 



( T c 



\10 4 K 



0.45 

I GHz/ 



-o.i 



15.1 



(— 

VGHz 



erg s" 



Hz- 



(15) 

This equation essentially weights the observed radio con- 
tinuum luminosity based on the expected thermal frac- 
tion at a given frequency. 

The most widely used radio SFR calibration, however, 
is the result of the tight, empirical FIR-radio correlation 
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TABLE 4 

Derived Parameters 



ID 




f 33 GHz 
JT 




Q NT t 




B*. 

min 


T c 






(%) 






(10 8 L o ) 


( M G) 


(10 4 K) 


Nucleus 


62± 6 


0.58±0.04 


0.74±0.05 


47.56± 2.14 


30. 3± 9.3 


0.42± 0.08 


Enuc. 


1 


93± 4 


0.10±0.01 


0.78±0.06 


1.08± 0.08 


10. 1± 3.0 


1.36± 0.31 


Enuc. 


2 


87± 5 


0.29±0.02 


0.82±0.06 


1.92± 0.11 


14.3± 4.1 


0.77± 0.27 


Enuc. 


3 


87± 7 


0.23±0.02 


0.79±0.06 


l.lOzt 0.08 


ll.Ozfc 3.2 


0.86± 0.31 


Enuc. 


1 


42*± 4 


0.20±0.01 


0.80±0.06 


1.32± 0.08 


11. 6± 3.4 


1.78± 0.59 


Enuc. 


5 


89± 8 


0.24±0.02 


0.77±0.05 


0.41± 0.05 


9.9± 2.9 


1.07± 0.44 


Enuc. 


6 


86± 5 


0.38±0.03 


0.83±0.06 


4.72± 0.25 


16. 0± 4.6 


1.04± 0.23 


Enuc. 


7 


85± 5 


0.23±0.02 


0.83±0.06 


4.86± 0.26 


14.5± 4.2 


1.14± 0.25 


Enuc. 


8 


88± 5 


0.18±0.01 


0.81±0.06 


2.66± 0.14 


13. 1± 3.8 


2.80± 0.54 


Enuc. 


9 


85± 5 


0.28±0.02 


0.83±0.06 


4.12± 0.21 


14.3± 4.1 


1.22± 0.26 


Note. 




^ The observed and nor 


i-thcrnial radio spectral indices measured between 1.4 



a ii( I 8.5 GHz. * Calculated using the t otal 1 .4 GHz flux densities and the revised minimum 
energy calculation of Beck & Krausc (2005). * Calculated using the observed 33 GHz flux 
density before subtracting out the ~56% contribution thought to arise from anomalous 
dust emission (see fl2.4|l ■ 



(|de Jong et al.lll985t iHelou et alJll985f >: 



Qm = l°g 



3.75 x 10 12 Li.4gh 2 



(16) 



While classically expressed using the integrated {IRAS- 
based) FIR luminosity, we use the total IR luminosity 
resulting in an average value of qir = 2.64 ±0.26 dex for 
a sample o f star-forming galaxies in the local universe 
(|Bellll206"3l) : this value is 0.30 dex larger than the aver- 
age value o f 2.34 ± 0.26 dex found when using the FIR 
luminosity ()Yun et al.ll200lD . Thus, one can express the 
1.4 GHz luminosity as a SFR by combining Equation 1161 
with Equation |4j 



/ SFRi.4GHz 

V M Q yr-i 



6.35 x 10" 



Li 



4GHz 



erg s _1 Hz" 



(17) 



It is worth stating here that this calibration is based 
on the FIR-radio correlation, which was established for 
globally integrated FIR and radio continuum properties 
of galaxies, and thus applies to galaxies as a whole. Most 
of the non-thermal radio emission from galaxies arises 
from electrons that have diffused away from their accel- 
eration sites in supernova remnants into quiescent regions 
void of current star formation. In fact, within individ- 
ual Hn region s, there may be little non-thermal emission 
present (e.g., iCondon fc Yinlll990[ ). Thus, locally near 
star-forming complexes, we expect the non-thermal emis- 
sion per unit star formation to be low given such a cali- 
bration as sug gested by studie s of the local FIR-radio cor- 
relation (e.g. JMurphv et alj l2006ai iHughes et all 120061: 
iTabatabaei et alj|2007j )~ 

4. COMPARISON OF SFR DIAGNOSTICS 

In the following section we compare SFRs estimated 
via the various diagnostics described in As previ- 
ously mentioned, we only include SFR diagnostics whose 
reliability is considered to be fairly insensitive to extinc- 
tion. SFRs estimated by free-free emission at 33 GHz are 
assumed to be the most accurate, and thus act as the ref- 
erence SFR for the comparison. This assumption is made 
as the thermal fraction at 33 GHz is found to be «87%, 
on average. The various SFR diagnostics are plotted in 
Figure [3] from left to right against the 33 GHz free-free, 



non-thermal, and total radio SFRs, respectively. The 
median ratio of the ordinate to abscissa (<R>), along 
with the residual dispersion after subtracting out an or- 
dinary least squares fit to the data (cr n t), are given in the 
top left corners of each panel. 

4.1. Non-thermal Radio Emission: 1.4 GHz Emission 
and the FIR-Radio Correlation 

In the top 3 panels of Figure G2 SFRs derived us- 
ing the observed 1.4 GHz flux densities (along with the 
FIR-radio correlation) from each star-forming region are 
plotted against the three 33 GHz derived estimates. Ex- 
cept for the starbursting nucleus, the 1.4 GHz SFRs ap- 
pear to significantly underestimate (i.e., by a factor of 
~5) the actual SFR as measured by the free-free emis- 
sion at 33 GHz. As stated at the end of §3.2[ this dis- 
crepant behavior is expected given that the FIR-radio 
correlation applies to galaxies as a whole. The agree- 
ment for the nucleus likely arises from the fact that 
the 33 GHz thermal fraction (i.e., ~62%) is significantly 
smaller than that of the extranuclear star-forming re- 
gions, being much more similar to what is found fo r glob- 
ally i ntegrated measurements of galaxies (e.g., ICondon! 
1992). Differences between the radio and optical prop- 
erties (e.g., 1.4 GHz to Ha ratios) of galaxy nuclei and 
Hll regions ha s been noted in the past by a number of 
authors fe.g.. [Israel fc Kennicuttl Il980t iHeckman et al.l 
119831 : Ivan der Hulst et all 11988^ These differences ap- 
pear to be due to a larger non-thermal fraction in the 
nucleus relative to what is estimated for the extranu- 
clear r egions as previously suggested by iKennicutt et all 
(1989), rather than simply higher optical extinctions. 

Comparing the 1.4 GHz SFRs with the 33 GHz SFRs 
calculated using the total 33 GHz flux densities (and 
Equation [T5)) . we again find that the 1.4 GHz SFR cal- 
ibrations are underestimating the SFR as measured by 
the total 33 GHz flux density among the extranuclear re- 
gions. This is similar to the discrepancy between the 
1.4 GHz and 33 GHz free-free SFRs, and most likely due 
to the fact that the thermal fractions are large at 33 GHz 
among these star-forming regions. On the other hand, 
the 1.4 GHz SFRs are nearly equal to those measured 
using the non-thermal 33 GHz flux densities among both 
the extranuclear regions and the nucleus. This is not too 
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Fig. 3. — A comparison of the various extinction independent 
SFR diagnostics described in [[3] (i.e., 1.4 GHz via the FIR-radio 
correlation, total IR, 24 fim, Ha + 24 /im — mix, and UV + IR — 
tot), plotted from left to right against the 33 GHz free-free, non- 
thermal, and total radio SFRs, respectively. To help guide the eye, 
the one-to-one line is shown as a dotted line in each panel. The 
points are also identified per column in the bottom row of panels by 
the extranuclear region number and an 'N' for the nucleus. In the 
upper left corner of each panel, the median ratio of the ordinate 
to abscissa (<R>), along with the residual dispersion about an 
ordinary least squares fit to each trend (ctri), are given. In the plots 
comparing the Ha+24 fim (mix) and 33 GHz SFRs, we additionally 
plot the nucleus using t he Ha +24 fim calibration for entire galaxies 
(Kcnnicutt ct al. 2009) as a cross. 



surprising given that the 1.4 GHz flux density is dom- 
inated by non-thermal emission, and suggests that our 
thermal decomposition is in fact fairly robust. The only 
region where there is a significant difference is extranu- 
clear region 4, however, the thermal/non-thermal decom- 
position here is complicated by the presence of an anoma- 
lous emission component. 

4.2. Total Infrared 

The second row of panels in Figure [3] compares the 
IR- based SFRs to those at 33 GHz. We find that the IR- 
based SFRs typically underestimate the SFRs derived 
from the free-free emission among the extranuclear re- 
gions by an average factor of ^2. There may also be 



a slight trend of increased discrepancy between the two 
SFR indicators with decreasing infrared surface bright- 
ness. This discrepancy among the extranuclear regions 
likely arises from the subtraction of a local background 
together with the application of a SFR calibration that 
assumes all UV photons are absorbed by dust at the 
source, which may not be true, since a significant frac- 
tion of non-ionizing photons can escape the star-forming 
complexes, powering the diffuse interstellar dust emis- 
sion. That is, the lower surface brightness extranuclear 
star-forming regions may become increasingly porous to 
UV photons. Thus, by subtracting a local background, 
which makes up a significant fraction of the cool/cold 
dust emission, the IR-based SFR calibration yields val- 
ues that are significantly lower than the free-free emission 
estimates. Similar results are seen when comparing the 
IR and radio SFRs using the total 33 GHz flux densities. 

On the other hand, the nucleus is a factor of ~2 
times larger than the estimate using the free-free or total 
33 GHz emission. One reason for the large discrepancy 
in the nuclear SFR could arise from a significant amount 
of dust being heated by stars that are not contributing 
significantly to the ionization of H atoms. Alternatively, 
this difference may arise from the capture of ionizing pho- 
tons by dust, which would decrease the free-free emission 
at 33 GHz and increase the IR emission. A detailed dis- 
cussion about such scenarios is given in £15.11 

The IR-based SFRs are generally larger (i.e., by a fac- 
tor of more than ~2, on average) than the non-thermal 
33 GHz SFRs for all star-forming regions. This is con- 
sistent with the above 1.4 GHz SFR comparison among 
the extranuclear regions which suggest that the SFR cal- 
ibration based on the non-thermal radio emission typi- 
cally underestimates the true SFR in such regions. This 
finding is also in agreement with studies of the resolved 
FIR-radio correlation within galaxies that have reported 
elevated FIR/radio ratios relative to the canonical cor- 
relation value for individual star- forming regions (e.g., 
[ Murphy et aT1l2006al; iHughes et all 120061: iMurphv et all 
120081 ). The physical explanation for this observation is 
that CR electrons propagate significantly further than 
dust-heating photons, thus the non-thermal radio con- 
tinuum emission associated with individual star-forming 
regions is much more spread out and diffuse than the 
dust emission - an effect which is not taken into account 
by the SFR recipe (see ! }53|) . 

4.3. Warm Dust: 24 [im Emission 

Focusing on only the warm dust emission, we compare 
24 ^m-derived SFRs against the 33 GHz suite of SFRs 
in the third row of panels in Figure [3] Unlike the total 
IR-based SFRs, there is much better agreement between 
the 24 /xm-derived and free-free emission SFRs among 
all extranuclear star-forming regions. There is also rela- 
tively good agreement among the 24 /xm-derived and to- 
tal 33 GHz SFRs among all extranuclear star-forming re- 
gions, while the 24 /xm-derived SFR for the nucleus is a 
factor of ~2 larger than each of the 33 GHz SFR esti- 
mates. Similar to the discrepancy between the IR-based 
and 33 GHz free-free SFRs for the nucleus, we may ex- 
pect the 24 /xm-derived and 33 GHz free- free SFRs to also 
disagree because of either the capturing of ionizing pho- 
tons by dust, or enhanced heating of very small grains 
(see $5.1$ . Such grains are stochastically heated and only 
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depend on the strength of the radiation field energy den- 
sity. Warm dust may also arise from active galactic nuclei 
(AGN) activity, however there is currently no strong evi- 
dence for AGN a ctivity in the nucleus of NGC 6946 (e.g., 
iTsai et al J [2006( 1 . And, like other comparisons with the 
non-thermal SFRs, the 24 /zm derived SFRs are a fac- 
tor of ~3 larger than the non-thermal 33 GHz SFRs, on 
average. 

4.4. Extinction Corrected Ha: _ffa + 24 /im 

The hybrid Ha + 24 /mi SFR diagnostic has been found 
to be quite reliable relative to SFRs measured using 
extinction-corrected Paa recombination li ne observations 
(e.g HST NICMOS narrowband imaging; ICalzetti et all 
2007). In the fourth row of panels, we plot this SFR di- 
agnostic against each of the 33 GHz SFRs, keeping in 
mind that we believe that the free-free emission likely 
gives the most reliable SFR. We find a similar behavior 
here as with the SFRs derived using warm dust emission. 
The Ha + 24 fj,m SFRs agree quite well with the 33 GHz 
free-free SFRs for the extranuclear star forming regions, 
however, it is slightly larger than (i.e., 60% above) the 
33 GHz thermal estimate in the nucleus. We find the 
same behavior in the comparison between the Ha +24 /im 
and total 33 GHz SFRs. The non-thermal 33 GHz SFRs 
again underestimate the SFRs among the extranuclear 
star-forming regions relative to the Ha + 24 /im SFRs by 
an average factor of ~4. 

4.5. Bolometric Estimates: UV + IR 

Finally, we compare our last extinction-insensitive SFR 
diagnostic (i.e., UV + IR SFRs) with the 33 GHz SFRs 
in the bottom panels of Figure [3] In general, we find 
that the UV + IR SFRs are in good agreement with 
the 33 GHz free-free SFRs among the extranuclear star- 
forming regions. The case is marginally improved when 
comparing the UV + IR SFRs to the total 33 GHz SFRs. 
Not surprisingly, we again find that the 33 GHz non- 
thermal SFRs are quite discrepant, begin a factor of 
nearly ~3 lower than the UV + IR SFRs. And, in all 3 
comparisons, we find that the nucleus UV + IR SFR is 
a factor of >2 larger than the 33 GHz SFRs. 

5. DISCUSSION 

In the previous section we have quantitatively com- 
pared a number of SFR diagnostics whose reliability are 
considered to be largely insensitive to interstellar extinc- 
tion for 10 star-forming regions within the nearby galaxy, 
NGC 6946. Using these results, we attempt to construct 
a self-consistent picture to explain the observed discrep- 
ancies between the SFR estimates for the nuclear and 
extranuclear regions by using free- free, dust, and non- 
thermal emission processes as diagnostics. Additionally, 
we discuss how combinations of these diagnostics can 
be used to infer other physical characteristics from star- 
forming regions, and star-forming galaxies as a whole. 

5.1. Anomalous Behavior of the Nuclear Starburst 

One of the most striking differences in the present com- 
parison among these star-forming regions is the behavior 
between the nucleus and the extranuclear star-forming 
regions. While the 1.4 GHz radio continuum emission 
and the FIR-radio correlation yield a SFR that agrees 



well with the 33 GHz free-free SFR for the nucleus, we 
find that all other diagnostics that rely on some sort 
of dust emission yield SFRs that are significantly larger 
than that when using the 33 GHz free-free emission. 

As previously stated, one way to explain the dust- 
based SFRs being larger than the 33 GHz free-free SFR 
estimates in the nucleus is through dust capture of ion- 
izing photons, which will in turn increase the dust emis- 
sion while decreasing the Ha and free-free emission per 
unit star formation. The extinction in NGC 6946's nu- 
clear starbur st has been measured to be Ay ~ 4.3 mag 
(8'.'5~280 p c:lEngelbracht et al.fi"9 961 and A v ~ 4.6 mag 
(2'/7~90 pc: lQuillen fc Yukitall200iF using Balmer decre- 
ments. Similar values have been estimated using the 
9.7 /im silicate a bsorption feature in the central kpc 
(Ay ~ 5.0 mag: iSmith et al.l 12007( 1. while much larger 
estimates have been quoted for t he central ~60 pc usin g 
total gas masses (Ay ~ 100 mag: lSchinnerer et alJ 2006). 
indicating a highly obscured, compact component. Al- 
though the extinction is high towards the central star- 
burst, the effect of dust capturing ionizing photons may 
not significantly affect our measurements over the entire 
~0.8kpc area covered by the GBT beam unless the star- 
burst is highly concentrated within the central ~60 pc. 

Using our 24 /im image at its native resolution, we find 
that ^25% of this warm dust emission arises within the 
central 200 pc. Assuming that 50 to 100% of this emis- 
sion from the central starburst region (i.e., 12.5 to 25% 
of the total 24 /im emission) is powered by ionizing pho- 
tons, which are thus not participating in the production 
of free-free emission, this does not account for the fac- 
tor of ~2 discrepancy between the SFR estimates as the 
24 /im SFR would still be larger than the 33 GHz free- free 
SFR by ^55 to 20%, respectively. This estimate relies 
on the relation betwee n the 24 /tm and free- fre e emission 
given by Equation 2 in iMurphy et al.l ((2006b) . It is also 
worth pointing out that the non-thermal (e.g., 1.4 GHz) 
emission should not be affected by this competition for 
ionizing photons and that these SFRs (i.e., the 1.4 GHz 
and 33 GHz non-thermal) agree with the 33 GHz free- 
free SFR. Thus, it seems that dust absorption of ionizing 
photons by itself cannot explain this discrepancy. 

Alternatively, it may be more likely that these rela- 
tively larger dust-based SFRs are the result of extra dust 
heating in the nucleus not associated with current star 
formation activity due to its large stellar density. In this 
case, the dust-based SFR will overestimate the true SFR 
for the nucleus of NGC 6946. Such a scenario seems plau- 
sible given the large stellar densities associated w ith the 
pseudobulge of NGC 6946 (|Kormendv et al.ll2010D which 
compl etely fills our beam; the radius of the bulge is 
»15 " dRegan fc Vog el 1995) and contains a stellar mini- 
bar (jElmegreen et al]ll9 98j. Furthermore, excess heat- 
ing of dust (small and large grains) may arise from an 
accumulation of non-ionizing stars in the nucleus as the 
result of an extended duration of star formation activ- 
ity associated with the starburst. Spectroscopic stud- 
ies of near-infrared recombination lines have found that 
star formation episodes in nuclei differ significantly from 
those in extranuclear Hn regions. Unlike extranuclear 
star-forming regions, galaxy nuclei typically show a spa- 
tial coincidence between the continuum and line emission 
suggesting that nuclear star formation activity is typi- 
cally ongoing, where as the star-forming episodes in ex- 
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Fig. 4. — The ratio of the 33 GHz free-free to 1.4 GHz (non-thermal) SFRs plotted against [/ ra d (left panel) and 1.4 GHz surface brightness 
(right panel). Ordinary least squares fits to the trends are shown. 



tranuclear regions are much younger, and more transient 
in nature (jDale et al.ll2004D . 

5.1.1. Galaxy Nuclei Behaving More Like Entire Galaxies 

Assuming that ionizing photons are not being captured 
by dust, it is interesting to point out that if we instead 
extinction co rrect the Ha line emis sion using the rela- 
tion given bv lKennicutt et all ((2009) for entire galaxies, 
which scales the 24 /im emission by a lower value (0.020 
versus 0.031 for Hn regions), we find good agreement be- 
tween the Ha + 24 ^m and 33 GHz free- free SFRs (crosses 
plotted in Figure Thus, it appears that the nucleus 
may behave more like an entire galaxy with respect to the 
non-thermal radio and Ha + 24/im SFR calibrations. As- 
suming then, that our UV + IR SFR calibration is only 
appropriate for extranuclear star-forming complexes, we 
can set our 33 GHz free- free SFR equal to Equation [5] to 
derive a new scaling to correct down the UV + IR SFR 
for the nucleus such that, 



SFRf " t c 
M Q yr" 1 



4.42 x 10 



-44 



L FUV + 0.43L IR 



erg s 



(18) 



The coefficient scaling the IR emission is actually more 
in line with a detailed study of UV and IR emission 
from global measurements of normal star-forming galax- 
ies, which report a value of ^0.46 (C. Hao, et al. 2011, 
in preparation). 

Admittedly, the fact that the FIR-radio correlation ap- 
pears to hold for the nucleus, and the non-thermal radio 
(i.e., 1.4 GHz and 33 GHz non-thermal) SFRs agree with 
the 33 GHz free-free SFR here, could be fortuitous. Ul- 
timately, it seems impossible to completely disentangle 
these two scenarios as there is no alternative way to ac- 
curately measure the ionizing photon rate. Thus, for such 
highly obscured regions, it appears difficult to accurately 
measure SFRs to better than a factor of ~2. 

5.2. Non-thermal Emission, SFRs, and CRs 

The fact that the prescriptions for converting the non- 
thermal emission into a SFR are doing a reasonable job 
for the nucleus relative to the extranuclear star-forming 
regions may be the result of the time dependence on CR 



electron diffusion. Work on the propagation of CR elec- 
trons in external galaxies has found a trend of decreasing 
CR electron propagation length with increasing radiation 
field energy density, as measured by the infrared surface 
brightness (jMurphv et all I2006U |2008[ ). The radiation 
field energy density in the nucleus, which appears to be in 
rough cquipartition with the magnetic field energy den- 
sity, is a factor of nearly ^25 times larger than that of 
the extranuclear star-forming regions, on average. This 
is consistent with a scenario of underestimated 1.4 GHz 
SFRs in the extranuclear star-forming regions relative to 
the nucleus. We can see if this scenario is valid by com- 
paring the cooling efficiencies and diffusion timescales for 
the nucleus and extranuclear star-forming regions. 

Let us assume that the propagation of CR elec- 
trons is described by a random walk with an energy- 
dcpendent diffusion coefficient De (e.g.. iGinzburg et all 
[1983), where D E = 5 x 10 28 cm 2 s" 1 for E < 1 GeV 
and 5 x 10 28 (£/GeV) 1/2 cm 2 s" 1 for E > 1 GeV (e.g., 
Jones et al J 1200 It IMoskalenko et al.l l2002TlMaurin etail 
20021) . Then, the time it takes to diffuse a distance Zdiff 



is given by Tdiff = 1% s /De- For 1.4 GHz emitting CR 
electrons to travel a distance of «0.8kpc (i.e., to cross 
the GBT beam at 33 GHz), it will take ~3.0 and 2.4 Myr 
in the nucleus and extranuclear regions, on average, re- 
spectively. The slight difference in Tdis arises from differ- 
ences in the typical energy of 1.4 GHz emitting electrons 
in the nucleus and extranuclear regions due to differences 
in the estimated magnetic field strengths. Accordingly, 
CR electrons in extranuclear star-forming regions will 
rapidly propagate away from their acceleration sites on 
timescales that are much longer than (i.e., a factor of ^5) 
estimates of their typical cooling times (~12 Myr). On 
the other hand, the estimated cooling time for 1.4 GHz 
emitting CR electrons in the nucleus is ^2.0 Myr, which 
is ~1.0 Myr shorter than the estimated escape time. 

Thus, the agreement between between the non-thermal 
and 33 GHz free-free SFRs for the nu cleus l ikely arises 
from being in the calorimeter limit of I Volkl (|1989h such 
that the continuously injected CR electrons, associated 
with the extended duration of star formation in the nu- 
cleus as discussed in H5.ll lose the bulk of their energy 
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before propagating outside of the region covered by our 
beam. It is also worth noting that the large gas densities 
in the nucleus will lead to an increase in hadronic interac- 
tions between CR nuclei and the interstellar gas, netting 
more synchrotron radiation per unit star formation due 
to an increase in pionic secondaries (e^). On the other 
hand, the transient nature of the star formation episodes 
in the young extranuclear star-forming complexes pro- 
vides enough time for CR electrons to propagate signif- 
icantly further (i.e., on the order of a ^kpc) than dust- 
heating and ionizing photons, resulting in an underesti- 
mate of the true SFR. 



5.2.1. Revised Non-Thermal SFR Recipes for Individual 
Star-Forming Regions 

It may be possible to correct the non-thermal SFR 
recipes by including a term that takes into account 
the time dependence of CR propagation. For example, 
through a detailed comparison between the non-thermal 
radio and infrared morphologies of nearby disk galax- 
ies, it has been shown that C/ ra( j can be used as a proxy 
for the average distance tr aveled by the CR electrons 
(jMurphv et al.ll2006U l2003). We therefore plot the ratio 
of the 33 GHz free-free to 1.4 GHz (non-thermal) SFRs 
against {/ ra d (left panel of Figure H]) . We apply an or- 
dinary least squares fit to derive a rough correction for 
the two SFR recipes that rely on the non-thermal radio 
emission (i.e., Equations fT4l and [17]) . Accordingly, we can 
re- write Equation [14] such that 



/sfr; 



NT-corr 



5.89lJ 5 5 x 10 



-32 



(-Y 

VGHz/ 

t^rad 



^erg cm 
and Equation [17] such that 



erg s _1 Hz -1 

-0.316±0.102 



(19) 



SFR? 



COIT 

4GHz 



M Q yr- 



5.64±|^ 3 x 10" 



32 



4GHz 



erg s 1 Hz 



ad 



erg cm 



-3 



-0.316±0.102 



(20) 



Uncertainties from the least squares fits are 
given. Using these revised expressions results 
in < SFRN t g - h 77SFRJ 3GHz > ~ 1.19 ± 0.36 and 
< SFR^ GHz /SFRT GHz > „ i.o7 ±0.33. 

Since U Ia d should roughly scale wi th Ub (e-g-, giv en the 
tightness of FIR-radio correlation; iCondonl 11992) . and 
Ub scales as the non-thermal radio surface brightness, 
it is likely more convenient to try to derive a correction 
which relies on the non-thermal radio emission alone. We 
therefore plot the ratio of the 33 GHz free-free to 1.4 GHz 
(non-thermal) SFRs against the 1.4 GHz surface bright- 
ness in the right panel of Figure 2] and perform a least 
squares fit. This results in corrections to the SFR recipes 
relying on the non-thermal radio continuum where Equa- 



tion 03 becomes 

/ CT?X> NT-corr \ 

f 2^ ) „ 7.33 +64 10 -36 

V M yr-i J - 6 ' 6 



_ NT / rNT 

f—Y ( » 

VGHz/ V WHz_1 



A.4GH 



erg s _1 cm -2 Hz -1 sr~ 4 



-0.374±0.057 



(21) 



and Equation [T71 becomes 

(exprjeorr \ 
„ ? ) ~ 7.01+g 1 , x 10- 35 
M Q yr-i ; - 6 - 3 



Lx. 



4GHz 



erg s 



Hz- 
h 



4GHz 



erg s cm 2 Hz 1 sr 



-0.374±0.057 



(22) 



These revised 

>-33GHz 

7SFR 



expressions result in 
■33GHz > ~ 1-15 ± 0.32 and 
< SFR™ r GHz /SFRj 3 GTJ , > - 0.98 ± 0.22. Accord- 



< SFR? T - corr 



ingly, these revised recipes may be more appropriate 
when estimating SFRs from individual star forming 
regions within galaxies that have had local background 
emission removed. 

5.3. Physical Parameters 

If we assume that Equation [6] provides accurate 
(extinction-corrected) Ha fluxes, and that the thermal 
radio continuum has been precisely quantified, we can 
compare these two quantities to obtain a rough idea of 
the electron temperature in each star-forming complex 
such that, 



To 



10 4 K 



3.68 x 10 23 



erg s 1 cm 2 Hz 



(23) 



f con 
J Ha 



erg s 1 cm 2 



(gHz) 



0.1 



-i 1.695 



(|Caplan fc Deharvend U986t iCondonl [1992T) . This re- 
lation assumes an ionized He to H density ratio of 
7V(He + )/iV(H+) - 0.08. The approximation for the 
temperature dependence of recombination line strengths 
are good to ~1% over electron densities and temper- 
atures ranging between 10 2 < n c < 10 4 cm" 3 and 
5000 < T e < 10000 K, respectively (jHummer fc Storey] 
Il987t ICondonl [1992). Under these assumptions, the scat- 
ter measured by plotting SFR m j x versus SFRjgQjjj, in 
Figure [3] is dominated by variations in the electron tem- 
peratures from each region. 

The derived electron temperatures (see Table range 
between 4000 < T < 14000 K, with a median of 
^11000 K, for all star-forming regions except for extranu- 
clear regions 4 and 8 where unrealistically high temper- 
atures are found, being 17000 and 28000 K, respectively. 
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- The same as Figure [3] except that SFRs have been 
using the photometry without local background sub- 



There are two possible reasons for such an outcome. One 
possibility is that the empirical Ha extinction correction 
could be underestimating the true Ha luminosity in these 
regions. This scenario seems plausible given that the ex- 
tinction correction is built upon an empirical correlation 
which has scatter, and different calibration coefficients 
for different luminosity ranges. 

Alternatively, the thermal fraction from these regions 
may be overestimated. We note that the 33 GHz flux 
density of extranuclear region 8 was boosted due to po- 
tential over-subtraction of background emission which 
may not have been adequate. For extranuclear region 
4, the situation may be more complex as the 33 GHz flux 
density from this source has been found to be a factor of 
^2 larger than what one would estimate by fitting the 
lower frequency (i.e., <10GHz) radio data. This excess 
emission is thought to arise from rapidly spinning ultra- 
small grains having a non-zero electric dipole moment. 
It is possible that the subtraction of this component was 
underestimated, however, the presence of such an excess 
indicates a high column of dust, which may suggest a 
possible underestimate in the Ha extinction correction. 
It is also worth noting that an over prediction of the ion- 



izing flux in extranuclear region 4 by the 33 GHz free- free 
emission (i.e., after subtracting out the excess emission) 
relative to the Ha + 24 /zm, 24 /im, and IR observations 
(e.g., see Figure U]) is opposite to the expectation if the 
excess emission was the result of optically-thick free-free 
emission arising from a compact H11 region. Presently, it 
is not possible to distinguish between these two scenarios. 

It is also worth pointing out the low derived electron 
temperature for the Nucleus. Similar to the discussion 
in §5.1, if we extinction correct the Ha using the rela- 
tion given by IKennicutt et a h (2009) for entire galaxies, 
the electron temperature for the nucleus appears more 
reasonable, being ~8000K. 

5.4. Applicability to Global Measurements of 
Star-Forming Galaxies 

While the current analysis has focused on individual 
star-forming complexes on ^0.8 kpc scales after removing 
a local diffuse background, it is instructive to investigate 
the significance of this subtraction in light of integrated 
measurements of star-forming galaxies for which a local 
subtraction is not possible. In Figure [5] we plot the same 
SFR comparisons as shown in Figure El except that the 
IR and radio SED fitting was done using the photome- 
try without subtracting out a local diffuse background 
component (see the Appendix for the associated fits and 
derived parameters). To briefly summarize the results 
from fitting the radio data in this manner, the observed 
radio spectral indices are steeper, being «0.5 on average, 
compared to «0.2 on average, when a local background 
was removed. The 33 GHz emission is still clearly dom- 
inated by free-free emission as the thermal fraction at 
33 GHz is only slightly smaller when including the local 
background emission, being 79%, on average. 

The trends found using these SFR estimates are nearly 
identical to those shown in Figure [3] The SFRs esti- 
mated using the non-thermal emission are still signifi- 
cantly underestimated relative to the 33 GHz free-free 
SFRs among the extranuclear regions, albeit by an av- 
erage factor of ^2 instead of ^4—5. This reduction in 
discrepancy is certainly due to the fact that the bulk of 
the non-thermal radio emission is diffuse. It then fol- 
lows that the corrections derived for the non-thermal 
SFR recipes in ^5.21 will only be applicable to discrete 
star-forming regions for which a local background has 
been removed. Refitting these trends using the quanti- 
ties derived before background subtraction leads to new 
coefficients for the SFR recipes in Equations [T9l EQl [2lJ 
and OH] of of 5.35 x 10" 31 , 5.11 x 1CT 31 , 3.31 x 10~ 33 , 
and 3.16 x 10 -33 , respectively, while the exponents scal- 
ing l^ad and Ii.4ghz become —0.209 and —0.267, respec- 
tively (see the Appendix for the re- written expressions). 
Interestingly, we again find good agreement between the 
non-thermal SFRs and the 33 GHz free-free SFR for the 
nucleus. 

Among the SFRs which rely on some form of dust 
emission, we again find a reasonable agreement with the 
33 GHz free-free SFRs for the extranuclear star-forming 
regions. One notable difference is that the IR-based 
SFRs appear to be in much better agreement with the 
free-free emission SFRs by not having subtracted a lo- 
cal background. We also find that the nuclear SFRs are 
larger than the 33 GHz free-free SFR for the nucleus by 
a factor of ~2.5, on average. So, while the sample pre- 
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sented here has targeted bright, active star- forming re- 
gions, and is by no means representative of the global 
emission from a normal star-forming galaxy, the ^0.8 kpc 
scales investigated here do average over a number of dis- 
crete star formation sites which likely have a range in 
ages. Given this, and how well these diagnostics seem to 
agree even when including local diffuse emission, it ap- 
pears that applying these diagnostics to entire systems 
which are vigorously forming star s is appropria t e. In 
fact, it has recently been shown bv lCalzetti et al.l (|2010l ) 
that the Ha + 24 /im calibration used here, and derived 
for Hn regions, works well for local IR-bright starbursts. 
Similarly, if we instead use the Ha + 24 fim calibratio n 
derived for entire galaxies (i.e., iKennicutt et al.l [2009n . 
which has been found to be appropriate for the 24 fim 
luminosities of each region (|Calzetti et al.l I2010D . when 
local backgrounds are not removed, these SFRs agree ex- 
tremely well with the 33 GHz free- free SFRs (i.e., median 
ratio near unity) over the entire luminosity range inves- 
tigated. Of course, applying such diagnostics to entire 
galaxies must be done with care. 

5.5. A Test for IMF Variations in High-z Galaxies 

Recent work on low surface brightness dwarf galaxies 
have uncovered discrepancies between SFR diagnostics 
which depend on the ionizing photon rate (i .e., extinction 
corre c ted Ha) and the UV continuum (e.g. JMeurer et all 
l2009t iLee et all I2009T I . One scenario to explain these 
discrepancies is a change in the upper end of the IMF 
of these systems since the ionizing photon rate is much 
more sensitive to changes in the high-mass end of the 
IMF than the UV continuum. Following the same logic, 
it may be possible to measure such discrepancies, and 
identify variations in the upper end of the IMF, in high- 
z objects using future high frequency {v > 10 GHz) ra- 
dio surveys, although it is worth noting that such vari- 
ations can be mimicked by age effects, stochastic sam- 
pling of the IMF, and ionizing photon escape. The radio 
may be particularly powerful since, as we have shown, at 
higher frequencies the radio emission from galaxies be- 
comes dominated by free-free emission. Surveys in the X- 
band (8 — 12 GHz) can probe the rest-frame 30 GHz emis- 
sion from star- forming galaxies at z ~ 2. In addition to 
measuring higher rest-frame frequencies with increasing 
redshift, the non-thermal component of a galaxy's radio 
continuum emission is expected to become increasingly 
suppressed with increasing redshift due to rapid cooling 
of CR electron s off of the cosmic microwave background 
(|Murphvl [2009) . Hints of such an effect have been ob- 
served (e .g., fiat radio spec tral indices of l ensed, high- z, 
galaxies; llvison et all 120101: llkarashi et all 120101 ) . This 
additional suppression of the non-thermal emission from 
galaxies could improve the situation by leaving only the 
thermal component detectable. 

In Figure £0 we present the same comparison of SFR 
diagnostics shown in the bottom 3 panels in Figure [5l 
except that top-heavy IMF calibrations have been used. 
We again run Starburst99 under the same assumptions 
as described in [J3j except that the slope of the upper 
mass function has been flattened from —2.3 to —1.5. We 
choose a top heavy slope of —1.5 as this has been sug- 
gested as a way to explain the reionization of the in- 
tergalactic m edium at z < 11 by star- forming galaxies 
(|Chary||2008[ ). In the associated output the ionizing flux 
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Fig. 6. — The same as the bottom 3 panels in Figure [5] except 
that a top-heavy Kroupa IMF was used, where the slope of the mass 
function is -1.5 for stellar masses ranging between 0.5 — 100 Mq. 

increases by a factor of «6 by moving to a top-heavy 
IMF, while the UV and IR emission only increases by 
a factor of «3. One of course expects that SFRs which 
depend only on the ionizing flux should be similarly inde- 
pendent of the IMF. However, based on these differences 
in the calibrations, SFRs proportional to the ionizing flux 
should be a factor of ~2 times larger than those which 
depend on the UV and/or IR luminosities for a top heavy 
IMF compared to a standard Kroupa IMF. In looking at 
Figure [6] this behavior is seen; the UV + IR SFRs yield 
a SFR which is nearly a factor of ^2 discrepant from the 
33 GHz thermal radio continuum estimate (i.e., a factor 
of ^2 larger since the IMF of NGC 6946 is more like a 
standard Kroupa IMF). Similarly, a systematic increase 
of ^40% is found when comparing the IR+UV SFRs to 
the SFRs estimated from the total 33 GHz radio contin- 
uum. 

We note that in this last comparison, all three quan- 
tities used to estimate the SFRs are as observed; local 
backgrounds were not subtracted, no extinction correc- 
tion was necessary, nor does it depend on an empirical 
calibration. Thus, by having high-quality UV, IR, and 
thermal radio continuum emission in hand, one should be 
able to quantify such systematic discrepancies between 
SFR diagnostics, and search for sources with potentially 
top heavy IMFs. 

6. CONCLUSIONS 

In the present analysis we have investigated the ap- 
plicability and reliability of a number of SFR diagnos- 
tics that are not hampered by extinction from dust. 
They include 1.4 GHz (non-thermal) radio, IR, 24 /im, 
Ha + 24 /im, and UV -I- IR observations. This was done 
under the assumption that our free-free emission mea- 
surements in the Ka-band (i.e., at 33 GHz) provide highly 
accurate estimates for the ionizing flux of massive star 
forming regions. Although we only focus on 10 star- 
forming regions in a single galaxy (i.e., NGC 6946), the 
nuclear star-forming region is classified as a starburst and 
these regions span nearly 2 orders of magnitude in IR lu- 
minosity Our main conclusions can be summarized as 
follows. 

1 . Excluding the total IR and non-thermal radio emis- 
sion, there is very good agreement between the 
various diagnostics and the SFRs derived from the 
33 GHz free-free emission among the extranuclear 
star-forming regions, independent of local back- 
ground subtractions. While the IR-based SFRs 
typically underestimate the SFRs derived from the 
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free- free emission by an average factor of ~2 among 
the extranuclear star-forming regions, there is sig- 
nificantly better agreement between these two diag- 
nostics when local backgrounds are not subtracted. 
These differences likely arise from the subtraction 
of a local background and the application of a SFR 
calibration that assumes all UV photons are ab- 
sorbed by dust at the source. 

2. Each of the diagnostics that rely on some sort of 
dust emission component yield SFRs for the nu- 
cleus that are larger than the 33 GHz free- free emis- 
sion SFR by an average factor of ~2. Given that 
the 1.4 GHz SFR agrees with the 33 GHz free-free 
emission SFR, and that only ^25% of the 24 /im 
emission within our beam arises from the central 
starburst, it seems most likely that these discrep- 
ancies are the result of excess IR heating arising 
from an accumulation of non-ionizing stars associ- 
ated and extended duration of star formation in the 
nucleus rather than the capturing of ionizing pho- 
tons by dust. However, conclusively distinguishing 
between these two scenarios is difficult. 

3. Typical SFR recipes that rely on non-thermal radio 
continuum emission (e.g., 1.4 GHz SFRs based on 
the FIR-radio correlation) underestimate the true 
SFRs by an average factor of ^4—5, or an aver- 
age factor of ^2 when local backgrounds are not 
removed, among the extranuclear star-forming re- 
gions. This is in contrast with the nucleus, for 
which the non-thermal emission estimates agree 
rather well with the 33 GHz free-free emission. 
We believe that the agreement between the non- 
thermal and 33 GHz free-free SFRs arises from CR 
electrons decaying within the nuclear starburst re- 
gion with negligible escape, mirroring its high opti- 
cal depth to ionizing photons. On the other hand, 
the transient nature of star formation episodes in 
the young extranuclear star-forming complexes al- 
lows for CR electrons to diffuse significantly further 
than dust heating photons, resulting in an under- 
estimate in the SFR using the non-thermal emis- 
sion. To properly measure the SFRs from individ- 
ual star-forming regions using non-thermal emis- 
sion requires a revised SFR recipe that attempts 
to incorporate the time-dependence of CR propa- 
gation. 

4. Using the total 33 GHz flux densities (i.e., not 
corrected for the presence of non-thermal emis- 
sion) yields SFRs that are in excellent agreement 
with those measured using only the 33 GHz free- 
free emission and combination of Ha + 24 /im data 
among the extranuclear star-forming complexes. 
This result arises from the fact that free-free emis- 
sion dominates the total flux density at 33 GHz 
(i.e., «87%, on average, among the star- forming re- 



gions investigated) and suggests that observations 
in the rest-frame Ka-band may provide accurate 
SFRs in galaxies without the need of ancillary radio 
data to account for the non-thermal component. 

5. Even without correcting our photometry for local 
diffuse background emission, which is prominent 
for observations that are most sensitive to non- 
thermal radio or cold dust emission, the 33 GHz 
flux densities are still dominated by free-free emis- 
sion with thermal fractions of «79%, on aver- 
age. Furthermore, comparisons between the var- 
ious SFR diagnostics remain generally unaffected, 
suggesting that these recipes are applicable to en- 
tire galaxies which are dominated by ongoing star 
formation activity. 

6. Using a combination of high frequency {v > 
10 GHz) radio data along with rest-frame far- 
infrared and UV data, one can compare radio ver- 
sus UV + IR SFR estimates which should become 
increasingly discrepant (i.e., a factor of ~2) as the 
IMF shifts towards one that is top heavy. This 
test makes use of the fact that radio continuum 
emission becomes dominated by thermal (free-free) 
emission at higher frequencies and, in addition, 
that the non-thermal emission of high-z galaxies 
should become increasingly suppressed due to in- 
creased IC losses of CR electrons off of the cos- 
mic microwave background. Thus, deep surveys at 
frequencies > 10 GHz (e.g., X-band; 8 - 10 GHz) 
should begin to detect the free-free emission from 
galaxies at z > 2 directly. Such tests to identify 
variations in the IMF of high-z systems should be 
possible using deep survey imaging from upcoming 
radio (EVLA, MeerKAT & SKA/NAA), mm/IR 
(ALMA/CCAT), and optical/NIR (JWST) facili- 
ties. 
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APPENDIX 

ESTIMATES OF LOCAL BACKGROUNDS AND THEIR EFFECT IN THE SED FITTING 

In Table Owe provide the values used for the local background subtraction at each wavelength for all star-forming 
regions. We also give the average fractional contribution of the background emission to the total emission per each 
wavelength and region at the end of each row and column, respectively. These values provide an estimate of the 
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TABLE 5 
Local Background Measurements 



Band 


Nuc. 


Enuc. 1 


Enuc. 2 


Enuc. 3 


Enuc. 4 


Enuc. 5 


Enuc. 6 


Enuc. 7 


Enuc. 8 


Enuc. 9 


/Jiff ( % ) 


Si. 4, GHz (mjy) 


14.51 


2.98 


2.49 


2.00 


1.61 


1.43 


5.38 


5.63 


4.11 


4.56 


51 


Si. s GHz (mjy) 


13.82 


2.79 


2.17 


1.78 


1.34 


1.01 


5.17 


5.10 


3.87 


4.14 


49 


Si. 7 GHz (mjy) 


12.37 


2.49 


1.96 


1.56 


1.16 


0.81 


4.66 


4.60 


3.34 


3.74 


47 


04.9 GHz (mjy) 


5. 14 


0.90 


0.71 


0.64 


0.65 


0.61 


1.90 


1.77 


1.49 


1.48 


33 


■So c i~i tj ( tn Tvi 
^8.5 GHz V illu J/ 


3.53 


0.80 


0.55 


0.55 


0.38 


0.26 


1.49 


1.23 


1.04 


1.13 


31 


/4850Mm) (mjy) 


42.0 


1.2 






0.0 




4.3 


2.7 


0.0 


3.6 


5 


/„ (350 Atm) (Jy) 


0.99 


0.21 


0.16 


0.14 


0.17 


0.09 


0.48 


0.40 


0.31 


0.35 


45 


/4250/xm) (Jy) 


2.71 


0.52 


0.38 


0.36 


0.43 


0.22 


1.26 


1.05 


0.83 


0.87 


55 


jv(160 Aim (Jy) 


6.48 


1.03 


0.74 


0.54 


0.63 


0.40 


2.81 


2.31 


1.66 


1.83 


38 


MlOOAtm) (Jy) 


5.90 


0.85 


0.54 


0.40 


0.39 


0.24 


2.38 


1.89 


1.32 


1.45 


30 


/470/zm) (Jy) 


3.12 


0.37 


0.22 


0.17 


0.13 


0.06 


1.23 


0.89 


0.57 


0.66 


21 


/„(24 M m) (mjy) 


169.1 


24.7 


14.4 


12.1 


11.3 


7.8 


63.2 


55.0 


37.2 


41.9 


15 


/„(8/im) (mjy) 


170.6 


32.4 


18.3 


17.2 


13.4 


13.2 


77.2 


68.3 


50.8 


54.8 


28 


/Ha (10~ 12 ergs -1 cm~ 2 ) 


0.34 


0.24 


0.16 


0.07 


0.07 


0.07 


0.32 


0.21 


0.13 


0.15 


16 


/ v (1528 A) (mjy) 


0.49 


0.93 


0.69 


0.31 


0.41 


0.36 


1.30 


0.84 


0.71 


0.69 


37 


/42271A) (mjy) 


1.74 


1.88 


1.67 


0.62 


0.94 


0.83 


2.86 


2.01 


1.82 


1.67 


40 


/iff (96) 


14 


49 


21 


29 


22 


39 


40 


32 


46 


37 





Note. — ' The average fractional contribution of the background emission at each waveband and region (i.e., the fraction of diffuse emission) arc 
given at the end of each row and column, respectively. 



TABLE 6 

Derived Parameters Calculated without Local Backgrounds Subtracted 



ID 




f 33 GHz 


a obst 


Q NT t 




mm 


T e 






(%) 






(10 8 L o ) 


( M G) 


(10 4 K) 


Nucleus 


53± 6 


0.62±0.04 


0.74±0.05 


50.10± 2.43 


31. 6± 9.7 


0.31± 0.07 


Enuc. 


1 


80± 5 


0.51±0.04 


0.81±0.06 


1.76± 0.11 


13. 8± 4.0 


0.72± 0.21 


Enuc. 


2 


83± 6 


0.43±0.03 


0.83±0.06 


2.31± 0.14 


15.9± 4.6 


0.65± 0.24 


Enuc. 


3 


79± 7 


0.48±0.04 


0.81±0.06 


1.45± 0.10 


13. 4± 3.9 


0.65± 0.25 


Enuc. 


1 


45*± 5 


0.36±0.02 


0.81±0.06 


1.62± 0.10 


13. 4± 3.9 


1.75± 0.61 


Enuc. 


5 


78± 8 


0.49±0.04 


0.80±0.06 


0.64± 0.06 


12. 1± 3.5 


0.73± 0.32 


Enuc. 


6 


79± 6 


0.56±0.04 


0.84±0.06 


6.21± 0.30 


18. 4± 5.3 


0.73± 0.19 


Enuc. 


7 


79± 6 


0.48±0.04 


0.84±0.06 


6.08± 0.31 


17.6± 5.0 


0.90± 0.22 


Enuc. 


8 


82± 6 


0.46±0.03 


0.83±0.06 


3.49± 0.18 


16. 0± 4.6 


2.14± 0.48 


Enuc. 


9 


79± 6 


0.50±0.04 


0.84±0.06 


5.05± 0.26 


17.0± 4.9 


0.96± 0.23 


Note 




^ The observed and nor 


-thermal radio spectral indices measured between 1.4 



a ii( I 8.5 GHz. * Calculated using the t otal 1 .4 GHz flux densities and the revised minimum 
energy calculation of Beck & Krausc (2005). * Calculated using the observed 33 GHz flux 
density before subtracting out the ~50% contribution thought to arise from anomalous 
dust emission (sec S|2.4t ■ 



typical amount of diffuse emission contributing to each wavelength and region. Not surprisingly, the low frequency 
(i.e., 1.4, 1.5, and 1.7 GHz) radio observations, which are dominated by non-thermal emission, are comprised of ^50% 
of diffuse emission, on average. The fractional contribution of the background to the total measurements at 850 /xm is 
anomalously low, which is most likely due to the insensitivity of these data. Comparing the average diffuse emission 
fraction among each star- forming complex, we see that the nucleus has the lowest contribution, being a factor of ^2.5 
smaller at 14%, on average. 

Additionally, we show the modeled SEDs (Figure [7]) and give the associated derived parameters (Tabled of each star 
forming region using the photometry before subtracting out a local diffuse component. The observed radio spectral 
indices are much steeper by not subtracting out a local diffuse component, being «0.5, on average, compared to ~0.2, 
on average, when the backgrounds are removed. This finding is consistent with the low frequency non-thermal emission 
being dominated by a diffuse component. Also note worthy is the fact that the thermal fraction at 33 GHz is still 
large, being «79%, on average. 

As described in ^5.41 the corrections derived for the non-thermal SFR recipes in H5.21 which attempt to take into 
account the time dependence of CR propagation, will only be applicable to discrete star-forming regions for which 
a local background has been removed. By refitting these trends using the photometry without subtracting local 
backgrounds leads to different revised SFR recipes such that Equations [T9] and [20] become 

/SFR NT ~ corl '\ / V \a NT / L NT \ / U \ -0.209±0.094 

V Moyr- 1 ) ~ 5 - 35 -" 89 X 10_31 (ghJ (ergs-i Hz~i J (erg cm" 3 J (1) 
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and 

( SFR 14GHz \ „ 5 n+ 55 x 1Q -31 f ^l-4GHz W ggd (2) 

V M yr- 1 J - 4 - 7 \erg s" 1 Hz" 1 J \erg cm" 3 / ' K ' 

respectively. Using these revised expressions results in < SFR^j H c ° rr /SFRj 3GHz > ~ 0.95 ± 0.27 and 
< SFR^°4q H7 / SFR^p. H7 > ~ 1.03 ± 0.29. Similarly, by using the photometry before local background subtraction, 



Equations [21] and [22] become 

/CT?-RNT-corr\ . . „»t / j NT \ / j \ -0.267±0.080 

\ 01+75 33/ V \ a ( L v \ ( il.4GHz \ 



V M Q yr-i ) ~ 3 - 31 - 32 X 10 VCHzJ Hz-i J Ug s^ 1 cm" 2 Hz"* sr-i J (3) 
and 

(QTTTJcorr \ / t \ / r \ -0.267±0.080 

^^1.4GHz \ 3 1 6 +72 x 10 -33 / ^1.4GHz W A.4GHz \ f ) 

M yr-i J 6 -- 3 U l^erg s^ 1 Hz"* J I, erg s" 1 cm" 2 Hz-* nr-* J ' 1 j 

respectively. These revised expressions result in < SFR^3^° rr /SFR3 3GHz , > ~ 0.97 ± 0.24 and 
< SFRf | GHz /SFRg 3GHz > ~ 1.03 ± 0.25. Accordingly, these revised recipes may be more appropriate when 
estimating SFRs from individual star forming regions within galaxies that have not had local backgrounds removed. 
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